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Instanton size distributions from calibrated cooling 
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Using an under- relaxed cooling algorithm we investigated the vacuum in the 2d 0(3) model and 4d pure gauge 
SU(2). We calibrated the amount of cooling performed to have si milar physical eff ect at different lattice spacings. 
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1. Introduction 

Topology is an important aspect of many lat- 
tice field theories, in particular QCD, where a 
non-zero value of the topological susceptibility is 
necessary to solve the U(l) problem and explain 
the mass of the rj meson. For QCD the topol- 
ogy is a consequence of the existence of instanton 
solutions of the theory, so simpler theories than 
QCD that share the property of instanton solu- 
tions have long been of interest. 

In this paper we present work carried out on 
two such simpler theories: 2d 0(3) and Ad pure 
gauge SU(2). In both of these theories we ex- 
tracted information about the underlying long- 
range structure of the theory by cooling, a pro- 
cess which has been widely used for this pur- 
pose. It was immediately apparent to us that, as 
we wished to discuss the underlying physics, we 
would need to compare results obtained at differ- 
ent lattice spacings and we would therefore need 
to calibrate the cooling we used to have the same 
physical effect across the range of a that we used; 
performing the same, arbitrary number of cool- 
ing sweeps at each value of a would have different 
physical effects at different couplings, rendering a 
comparison of the results meaningless — different 
cooling levels would mean the gross structure of 
the configuration would be different, and it is on 
this structure that we wish to make calculations 
and compare results. 



2. The 0(3) Model 

The 2d 0(3) model shares many characteristics 
with QCD: it is asymptotically free, with dynamic 
mass generation and instanton solutions due in 
this case to the homotopy class of windings from 
S2 — ► S^- 
Using a standard stereographic projection the 
single instanton solution is given by: 
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with both p and r complex, and <j>i the ith com- 
ponent of <p. The action of this continuum held 
configuration is Sj = Air/g 2 . (It should be noted 
that on a torus there is no single-instanton solu- 
tion p|.) In general multi- instanton solutions are 
given by analytic functions of u, ie c\lo = 0. 

Unfortunately, the lattice version of the 0(3) 
model has problems, as the theory is dominated 
by short-range fluctuations g], and formulating 
the theory on a lattice imposes a minimum size 
to these fluctuations, and so the important contri- 
bution from those smaller than a is lost. Further- 
more, additional unphysical contributions arise 
from lattice artefacts of size 0(a). So on the 
one hand results obtained from cooling would not 
have the contribution from the physical ultra- 
violet fluctuations, but on the other hand they 
would not be corrupted by the unphysical arte- 
facts, and measurements of objects of moderate 
size, a <C p <C L, should be reliable. On these 
grounds we decided to investigate the size distri- 
bution of instantons. 



Table 1 

Details of the cooling used for 0(3). Data calcu- 
lated on a 64 2 lattice. 
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C/3 £/2 2£/3 



0.80 
0.84 
1.00 



0.261(1) 
0.318(1) 
0.551(1) 



28 
17 
5 



79 
45 
13 



220 
106 

24 



We used an under-relaxed cooling given by: 



E, 



'->X+fJ, 



c-„) (2) 



with 4>' x normalised to unit length. This was used 
for two reasons: firstly, as a deterministic update 
it is much faster to compute, and secondly: the 
parameter a allows us to control the severity of 
the cooling, though a trade-off must be made be- 
tween the gentler cooling associated with larger 
values of a and the greater number of sweeps 
then needed to reach a given physical state. We 
eventually decided that a = 2 was a good com- 
promise between too many sweeps and too harsh 
a cool. O(a) effects in the lattice 0(3) action 
mean that lattice instantons are unstable under 
cooling and will be shrunk and eventually anni- 
hilated by prolonged cooling. We decided to cal- 
ibrate our cooling for the 0(3) model by mea- 
suring the number of cooling sweeps at different 
values of a required to annihilate a configuration 
generated from a discretised form of eq. |l|. We 
looked at three levels of cooling at each of three 
couplings: g 2 = 1.00,0.84,0.80 corresponding to 
mass gaps of ma = 0.551(1), 0.381(1), 0.261(1), 
the levels being those required to remove objects 
of size one-third, one half and two-thirds the cor- 
relation length £ = 1/m at each value of g 2 . The 
number of sweeps needed are given in table H. We 
found that S/SiV calculated on the cooled con- 
figurations was consistent across g 2 for each of the 
cooling levels we used, and took this as evidence 
for correct calibration of our cooling. 



3. Pure gauge SU{2) 

For SU(2) we again used an under-relaxed cool- 
ing, this time given by: 



U 'x.u = aU x^ + £*,/* 
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with Ui ,, normalised to lie in SU(2) and Yi x .„. the 



X,fl 



J x,i± 



sum of the 'staples' around U Xyll . 

When we came to extend the techniques we de- 
veloped for 0(3) to the case of Ad SU(2) we ran 
into the problem that the lattice SU(2) instan- 
tons are much more stable under cooling than 
their 0(3) counterparts — so stable in fact that 
had we naively taken the same criterion as ear- 
lier, we would have cooled the configurations be- 
yond any region where we would wish to study 
the vacuum, and in certain cases removed all the 
physics present in the uncooled configuration. In 
light of this we changed our approach: as we wish 
to consider the vacuum as quantum fluctuations 
around classical solutions, it makes sense, taking 
into account the stability of the lattice instantons, 
to calibrate the cooling by first constructing clas- 
sical solutions, introducing quantum fluctuations 
by a few Monte Carlo updates and then cooling 
until the original configuration is recovered. This 
procedure is performed a number of times and the 
average number of sweeps needed to recover the 
initial state is taken as the calibration. Prelim- 
inary studies of this on 16 4 lattices at (3 = 2.4 
and 24 4 at (3 = 2.5 have proved promising, and 
the data imply that 61 sweeps at a = 2 is suf- 
ficient to remove the quantum fluctuations with- 
out overly disturbing the underlying structure on 
a 16 4 lattice at j3 = 2.4. 

4. Instanton Size Distributions 

In order to calculate the sizes of instantons we 
looked for connected regions around local maxima 
in the action density for which the action den- 
sity was not less than half the value at the maxi- 
mum under consideration; 2-dimensional regions 
in 0(3) and 4-dimensional for SU{2). We then 
took the appropriate root of this volume to ob- 
tain a size parameter. Predictions of the instan- 
ton size distribution in both 0(3) and SU(2) can 
be derived from the dilute instanton gas model 
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predicts LTH-331, to appear in Phys. Rev. D. 

1 ^ l (4) 



SiV dp p 

whereas our data, presented in H and sum- 
marised in figure pi indicate a much stronger UV 
dependence and a distribution ~ p~ 3 . For SU(2), 
the dilute gas predicts an infra-red divergence and 

It is unclear whether our preliminary data, shown 
in figure [ll are in agreement with this. Certainly 
they indicate no ultra-violet divergence, and it 
appears that the very smallest instantons are all 
but absent. 

For 0(3), while the finite lattice spacing re- 
duces any signal below p ~ l-5£, we see a rapid 
decrease at larger p. For SU(2) this decrease is 
less rapid. The physical volumes here are approx- 
imately equal, with L ~ 16£, so we take this dif- 
ference in large p behaviour to be a physical effect. 

For 0(3) we also looked at the separation of 
pairs of instantons and anti-instantons, and found 
evidence for an interaction: the closest separa- 
tion of unlike pairs of objects was significantly 
smaller than that of like pairs. At each value of 
g 2 we found that unlike pairs had a closest separa- 
tion only 70% that of like pairs. This interaction, 
and the form of the distribution we found show 
no support for the dilute gas model for 0(3) in- 
stantons on the lattice. The corresponding SU(2) 
measurements are in progress. 

5. Conclusions 

We have given an argument in favour of cal- 
ibrating cooling to have a given physical effect, 
and altering the number of sweeps as the lattice 
spacing is altered. We have also presented data 
calculated using calibrated cooling and shown 
how it is consistent across a range of couplings. 

REFERENCES 

1. Solitons and Instantons, R. Rajaraman, 
North Holland (1982). 

2. M. Luscher, Nuc. Phys. B200, 61 (1982). 



0.020 



0.015 — 



□ ma = 0. 551(1) 

O ma = 0. 318(1) 

0.010 I— O ma = 0. 261(1) 



T3 \- o ma = 0. 318(1) 

> 



T3 



0.005 — 



0,000 



1 .0 M I I I I I I I I I I I I I I I I I I I I I 

0.9 E- 

0.8 — m 0+ =0.94(8) 

0.7 ^V=16 4 

0.6 — 



■§■ 0.5 - 
> = 

^ 0.4 E- 



0.3 — — 

; m : 

0.2 — m m ~ 

E ° m m E 

0,1 ^- mi 

Q Q Mlllll|jlllllllllllllllllllllllllllllllllllllllllllllllllir 

1 23456789 10 11 12 
P 



Figure 1. The distributions of instanton size 
calculated for 0(3) from 1000 configurations for 
each value of a on a 64 2 lattice (upper plot) , and 
for SU(2) from 1000 configurations at (3 = 2.4 
on a 16 4 lattice (lower plot). Instanton sizes are 
given in units of the correlation length, £ = 1/m; 
for 517(2), m = m +. We define iV = S/S* 7 



